Quantum dynamics in atomic-fountain experiments for 
measuring the electric dipole moment of the electron with improved sensitivity 
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An improved measurement of the electron electric dipole moment (EDM) appears feasible using 
ground-state alkali atoms in an atomic fountain in which a strong electric field, which couples to 
a conceivable electron dipole moment (EDM), is applied perpendicular to the fountain axis. In a 
practical fountain, the ratio of the atomic tensor Stark shift to the Zeeman shift is a factor fi ~ 100. 
We expand the complete time evolution operator in inverse powers of this ratio; complete results are 
presented for atoms of total spin F — 3, 4, and 5. For a specific set of entangled hyperfine sublevels 
(coherent states), potential systematic errors enter only as even powers of making the expansion 
rapidly convergent. The remaining EDM mimicking effects are further suppressed in a proposed 
double-differential setup, where the final state is interrogated in a differential laser configuration, 
and the direction of the strong electric field also is inverted. Estimates of the signal available at 
existing accelerator facilities indicate that the proposed apparatus offers the potential for a drastic 
improvement in EDM limits over existing measurements, and for constraining the parameter space 
of supersymmetric (SUSY) extensions of the Standard Model. 
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I. INTRODUCTION AND OVERVIEW 
A. Significance of the electron dipole moment 

A permanent electric dipole moment of the electron or 
of any other fundamental particle, or of an atom in an 
eigenstate of angular momentum, is possible only if the 
symmetries of both parity (V) and time reversal (T) are 
violated [1,2]. By the CVT theorem [3] , T violation is the 
equivalent oiCV violation, where C is charge conjugation. 
CP violation has been observed, but only in the quark 
sector and only in the decay of the neutral K and B 
mesons. The Cabbibo-Kobayashi-Maskawa (CKM) mass 
mixing matrix in the Standard Model is consistent with 
these observations, but does not produce a large enough 
CP violating effect to account for the excess of matter 
over anti-matter in the universe [4-6] . 

Extensions of the Standard Model generically con- 
tain new massive particles and new sources of CP vio- 
lation that give rise to an electric dipole moment (EDM) 
for the electron. In field theories, the interaction of 
an electron with the electromagnetic field is described 
by the field-theoretical interaction Hamiltonian Hj(t) = 
j d 3 r fi(t, r), where the interaction Hamiltonian density 
ri = rl{t,r) is (in international mksA (SI) units) 

H = ec^-y^ipA,, -> eci)Ti i (q)i)A li . (1.1) 

Here 7 M is a Dirac matrix, e = — |e| is the electron charge, 
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c is the speed of light, and we assume that the fermion 
field operator ip is normalized so that \ip\ 2 has physical di- 
mension of inverse volume. The second form in Eq. (1.1) 
includes radiative corrections, where the vertex function 
is denoted as T M . For an electron on the mass shell, the 
vertex function in an arbitrary field theory can be ex- 
pressed in terms of Lorentz and CPT invariant terms as 

T"(q) = F 1 (q 2 )^ + F 2 (q 2 )-^cr^ q „ 

2m c c 

+ W)^(7 A W-2 TOo c 7 V) , 

where a^ v = i [ 7 M , 7 I/ ]/2, and a is the contribution of the 
anapolc moment [7]. (The form factors are the F\ Dirac, 
the F2 Pauli, and the Fa anapole form factor; the CP-odd 
F3 term which leads to the EDM does not carry a spe- 
cial name.) The first and second terms in Eq. (1.2) con- 
serve C, P, and T separately. The weak interaction gives 
rise to the last, anapole term, which conserves T and CP, 
but violates both P and C separately; it is not considered 
in the following. The third term involves the form fac- 
tor F3 and is odd under P and T, but conserves C In 
the non-relativistic limit, the electron EDM gives rise to 
the following effective Hamiltonian for the interaction of 
the spin of a free electron with an electric field, 

h E UM = (-do) 3 ■ E , d e = -F 3 (0) P^- . (1.3) 

2m e c 

In the effective Hamiltonian for a bound electron, one 
has to replace 3 — > RF where F = (F x , F y , F z ) is the 
total angular momentum of the atom (including the nu- 
clear spin) and R is the enhancement factor [sec Eq. (1.6) 
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below]. We recall that |e| = — e is the elementary charge 
unit. The Pauli spin matrices t? are dimensionless. The 
electric field E is multiplied by the elementary charge 
and the electron Compton wavelength h/(m c c) to yield 
an interaction energy. In the absence of a special struc- 
ture that would otherwise constrain the form factor F3 to 
be small, field theories which extend the Standard Model 
generally contain a significant EDM. The special struc- 
ture in the Standard Model is that CP violation occurs 
only in a single phase in the quark mass mixing matrix. 
The leading Standard Model contribution to a lepton 
EDM is at the level of four loops [8]. Because the re- 
sulting electron EDM is far too small to be observed by 
any proposed experiment, there are in practice no Stan- 
dard Model effects to account for, so mere observation 
of an electron EDM is direct evidence of a new and non- 
Standard-Model source of CV violation [9] . 

The sensitivity of EDM measurements to new phenom- 
ena is far-reaching. In some recent calculations, an elec- 
tron EDM arises through a mechanism that produces the 
neutrino mass [10] or is sensitive to physics at energy 
scales that exceed 10 s GeV [11], or is sensitive to dark 
matter [12], or that is responsible for baryogenesis [13]. 
The present limit on the electron EDM already presses 
supersymmetry (SUSY), especially when that limit is 
combined with those on the neutron EDM and those 
on the EDMs of diamagnetic atoms. Present limits on 
the electron EDM [14, 15] are lower by a factor of 100 
than EDMs predicted by some SUSY models [16-21] with 
super-partner masses of 100 GeV and CV violating phases 
of order unity. Therefore, these SUSY models could be 
excluded. The present EDM limits are also not in com- 
plete agreement with models with one-TeV super-partner 
masses [22]. 

In this paper, we examine in detail a proposal to signif- 
icantly lower experimental limits on the electron EDM, 
by roughly two orders of magnitude in comparison to 
present limits [14, 15]. A nonvanishing EDM on this 
level would imply the existence of new physics beyond the 
Standard Model, or, alternatively, imply that currently 
favored extensions of the standard model need to be sig- 
nificantly altered. Within SUSY, an unexpectedly small 
EDM could be realized if one assumes that CV violating 
phases are much smaller than currently assumed, or that 
the masses of superpartners are much larger than cur- 
rently expected. Quite sophisticated models have been 
proposed with this notion in mind: E.g., in so-called 
split Supersymmetry [23-25], one assumes that only the 
masses of those superpartners that contribute to EDMs 
most significantly, are larger than expected, whereas the 
masses of other SUSY partners remain in the expected 
range. However, in general, SUSY offers no special rea- 
son for small CV violating phases or any good reason for 
most of the superpartner masses to remain small [26] . If 
more accurate experimental results still turn out to be 
compatible with a zero EDM, one will begin to exhaust 
some of the simpler remedies and push the theory towards 
constructions inconsistent with the original motivations 
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B. Experimental idea and overview 

Our proposed measurement scheme is based on the ob- 
servation that the interaction Hamiltonian (1.3) is pro- 
portional to <7 • E [or F ■ E for a bound electron, see 
Eq. (1.6) below]. The influence of an EDM of the elec- 
tron is largest when an atom is put in an intense, static 
electric field to evolve a long time, and an atomic foun- 
tain apparatus (the details of which will be explained be- 
low) is an essentially undisturbed environment in which 
this can be done. The inclusion of an atomic fountain in 
EDM experiments has not been discussed in the litera- 
ture to the best of our knowledge. However, it is a subtle 
matter to select an actual observable that is sensitive to 
an EDM, takes advantage of the properties of an atomic 
fountain, cancels many systematic effects, and can be re- 
alized using lasers alone. Our proposed scheme consists 
of four steps. 

[Step 1.] We start with heavy alkali atoms of half- 
integer nuclear spin and prepare an initial state that is a 
coherent superposition of states within the upper hyper- 
fine manifold of the electronic S'1/2 ground state. 

[Step 2.] In the fountain, the state of the atom evolves 
from the initial one under the influence of a strong elec- 
tric field but only weak static magnetic fields. The full 
Hamiltonian generating the dynamics is given below in 
Eq. (1.12). The quantum dynamics of the atom are the 
main concern of the current theoretical investigation. Be- 
cause of the analytic simple structure of the Hamiltonian 
given in Eq. (1.12), which involves a tensor Stark term 
and an EDM term, the dynamics can be described semi- 
analytically using time-ordered perturbation theory. 

[Step 3.] The analysis of the time-evolved state pro- 
ceeds in a region where the atom is irradiated by laser 
light tuned to the S 1/2 \F,M) -> Py\F' = f + 1/2, M') 
transition (with j' = 1/2 or 3/2 and F = 7 + 1/2), 
roughly as follows (further details are discussed in the 
main body of the article): After many cycles of absorp- 
tion and spontaneous emission, any sublevcl S\/2\F, 
either ends in the dark state of the upper hypcrfine level 
(F = I + 1/2) with probability pm or ends in the lower 
hyperfine level (F = I — 1/2) with probability 1—pM- If 
an atom ends in the dark state, the atom "survives" in 
the upper hyperfine level. The probability that the ac- 
tual time-evolved state survives is given as the sum over 
the "survival probabilities" of its separate components. 
This probability depends on the electron EDM and de- 
fines the observable S(6) [see Eq. (2.26) below], where 9 
is the rotation angle (about the z axis in Fig. 1 below) of 
the propagation direction of the linearly polarized anal- 
ysis laser with regard to the coordinate system of the 
atomic fountain apparatus (see Fig. 1 below for the x, 
y and z axes). We here assume that the atomic state 
is initially prepared with respect to a quantization axis 
adjusted to be essentially parallel to the direction of the 
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applied electric field (which is the z axis in Fig. 1 below) . 

[Step 4.] As a separate step, the survival probabil- 
ity signal S(9) is measured by fluorescing the Si/ 2 ,F = 
1+1/2 to P^/2,F' = 1 + 3/2 cycling transition, then 
pumping all atoms from the lower S1/2, F = I — 1/2 hy- 
perfine level also into the S\/2,F = 1 + 1/2 level and 
fluorescing the cycling transition again. The ratio of the 
number of photons scattered over a fixed time while the 
cycling transition was fluoresced then equals the desired 
probability. Throughout this paper, we use the phrase 
"optical pumping" when the purpose of hitting atoms 
with a laser is to alter the atomic state, and "fluoresc- 
ing" when the purpose is to scatter photons from the 
laser that are to be counted. 

We propose to do all of this in a double-differential 
setup, where the interrogation laser in step 3 is tilted 
by angles 9 and —9, and the difference of the two sur- 
vival probabilities P(9) = S(8) — S(—6) is taken [see also 
Eq. (2.28) below]. That difference is measured for oppo- 
site signs of the applied electric field and the difference is 
taken again. This leads to the observable P°(9) (the su- 
perscript standing for "odd") that appears in Eq. (3.42) 
below. This double-differential setup is key to our pro- 
posal and eliminates many systematic effects. 

A schematic diagram for the resulting experiment to 
measure the electron electric dipolc moment using an 
atomic fountain is laid out in Fig. f . Atoms are collected 
and cooled in a magneto-optical trap, a fully coherent ini- 
tial state is prepared (so that the use of a density matrix 
becomes unnecessary), and the atoms are launched ver- 
tically. The atoms enter a region shielded against static 
magnetic fields. Within the atomic fountain, the atoms 
then rise into and fall out of a large electric field. If an 
electron EDM exists, the initial quantum state (step 1) 
rotates by a small angle about the electric field axis, while 
in the electric field (step 2). While still in the magnet- 
ically shielded region, but outside the electric field, the 
atoms are analyzed by optical pumping (state analysis, 
step 3), which has the effect of encoding the small angle of 
rotation into a relative shift in the population of the up- 
per and lower hyperfine levels of the ground state. Those 
populations don't change if stray magnetic fields are sub- 
sequently applied, so atoms can be allowed to fall out of 
the magnetic shield before the populations are measured 
by optical pumping and counting of the scattered pho- 
tons (fluorescence, step 4). 

The system of optical pumping for state preparation 
uses the property that an atom illuminated with a laser 
tuned to an allowed transition with A_F = has a unique, 
known "dark" state from which excitation by the laser is 
forbidden. The dark state is a continuous function of the 
cllipticity of the laser polarization; for a linearly polarized 
laser with its axis of polarization parallel to the z axis, 
the dark states is |P, 0), and for a laser propagating in 
the z direction that is circularly polarized with positive 
helicity, the dark state is | F, i* 1 ) . For conceptual simplic- 
ity, we mainly study initial states which are the coherent 
superpositions (\F,M) ± \F,-M))/V2. 



coordinate 
system 



electric field 
plates 



state preparation 




magnetic 
shield 



state analysis 



(i - - - ^> fluorescence 



magneto-optical 
trap 



FIG. 1: Schematic diagram of an experiment to measure an 
electron EDM using an atomic fountain. In the diagram the 
width of the parabolic trajectory has been greatly exagger- 
ated. In first approximation, we assume that atoms rise and 
fall strictly along the y axis, and that the applied electric field 
is strictly parallel to the z axis. Atoms are exposed to lasers 
at the points on the trajectory marked by solid circles. 



The system of optical pumping for state analysis 
(step 3) is similar to that of state preparation (step 1). 
An atom in any state of the upper, F = 1 + 1/2 hyperfine 
level that is illuminated with a laser tuned to an allowed 
transition with AF = will either end in the dark state 
of the laser or in some state of the lower hyperfine level. 
In the simple cases of a laser that is polarized parallel 
to the z axis, or that propagates parallel to z and is 
circularly polarized, we define the probability pm that 
an atom originally in the state | F, ends in the dark 
state after multiple cycles of excitation to the Py state 
and spontaneous emission. For 7r polarization, the dark 
state is \F, 0), because (M = 0)-(M = 0) transitions 
are only allowed if the total angular momentum changes. 
The probabilities pu depend on the branching ratios for 
the various spontaneous transitions and can be found as 
simple fractions by solving the 2F + 1 equations for the 
evolution of the states of a given hyperfine level; the es- 
sential physics for the probabilities is illustrated for one 
particular transition in Fig. 2. In the calculation, the 
probability that an atom of nuclear spin I in a P-state 
with quantum numbers j',F',M' undergoes a radiative 
transition to an S1/2 state with quantum numbers j, F, M 
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FIG. 2: An atom, initially in the state S1/2I4, M\, under- 
goes a large number of cycles of absorption of laser photons, 
followed by spontaneous emission. We define the probabil- 
ity pm to be the probability that an atom, illuminated with 
a 7r-polarized laser tuned to the transition £1/2, F = 4 to 
P1/2, F — 4, ends up in the dark Si/a |4, 0) state instead of 
ending in some state in the S1/2, F = 3 hyperfine level, af- 
ter many cycles of absorption and spontaneous emission. In 
the upper half of the figure, the arrows show the pattern for 
the first spontaneous emission and re-excitation for an atom 
initially in state S1/2 4, — 2), and the histograms show how 
the probability of being in any given level evolves after 0, 1, 
2, or 3 spontaneous emissions. The probability of ending in 
the dark Si/ 2 14, 0) state involves processes like S1/2 4, —2) — > 



Pi 



/2| 



-2) -> S 1/2 |4,-l) ->■ Pi 



/2\% -1) -> 5*1/2 14, 0) and 
similar. Values for pm are in the graph in the lower half of 
the figure, and correspond to those given in the second column 
of Table I. 



is proportional to 



(2j' + l) 2 (2F + 1) 



J 
m 



F 
•1 



F' 

J 



(1.4) 



where standard notation is used for the 3j and 6j sym- 
bols. One then sets up and solves the rate equations. In 
the rate equations, one may ignore the small frequency 
difference between the upper and lower hyperfine mani- 
folds of the ground state Si/ 2 \l ± 1/2, M) indicated in 
Fig. 2, because the frequency difference to the upper state 
Pj> is much larger than the hyperfine splitting. The sets 
of values p that are achievable for the case F — 4 are 
listed in Tables I and II; these turn out to be rational 
numbers. Values for F = 3 and F = 5 are found in 
Appendix B of Ref . [27]. 

For the cases of linear and circular polarization, the 



TABLE I: For / — 7/2, the probability pm that an atom, ini- 
tially in the state Si/2|4, M), remains in the F — 4 hyperfine 
level after being pumped with light tuned to the transition 
S1/2, F — 4 — > Pji , F' that is linearly polarized with the axis 
of polarization parallel to z. 
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M = 4 
M = 3 
M = 2 
M = 1 
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M = -1 
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990 
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4700 
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4700 
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3788344 
11524649 

1205694 
11524649 
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11524649 

151263 



J 7' i2 

2585 
1092 
2585 
112 
517 
7(1 
517 

112 
517 

1092 
2585 

1792 
258.-) 
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839216 
4101383 

133868 
4101383 

15680 
4101383 

2450 
4101383 

15680 
4101383 

133868 
4101383 

839216 
4101383 

1 



TABLE II: For / = 7/2 the probability p M that an atom 
initially in the state Si/ 2 14, 4) remains in the F = 4 hyperfine 
level after being pumped with light tuned to the transition 
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Pji, F', that propagates in the z direction, 



and is circularly polarized with positive helicity. 
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46289600 
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'182:', 


529525457747 



probability that an atom, emerging from the atomic foun- 
tain in a general superposition of states ^2 -^m \ F, M) 
for amplitudes Am, ends in the dark state is simply 
^2 \Am\ Pm- Such a probability encodes information 
about the original atomic state in question and can be 
easily measured using fluorescent detection. Final read- 
out and normalization to the incoming flux of atoms pro- 
ceeds as follows (step 4): A cycling transition from the 
S1/2 ground state with spin F = 1+1/2 to the P 3 / 2 state 
with F = / + 3/2is rung and the photons counted; then 
a laser tuned to the transition from the lower hyperfine 
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level of the ground state with F = I — 1/2 to a P state 
with F = 1+1/2 puts all atoms into the 5-state hypcrfinc 
level with F = I + 1/2, and then the cycling transition 
S1/2 — > P3/2 is fluoresced and the photons counted again; 
the ratio of the counts gives the probability. 

An alkali atom of large nuclear charge would be indi- 
cated for our atomic fountain because rclativistic effects 
enhance the EDM of such an atom compared to that of 
the free-electron by a large multiplicative factor R (the 
"enhancement factor"). The theory of the enhancement 
factor is well-established, and has been used to set the 
EDM limits in Ref. [28]. Computed values of the en- 
hancement factors for Rb, Cs, and Fr are presented in 
Table III. None of these computed factors have varied 
by more than 20% from the earliest factors computed 
in 1966 (Ref. [29]), and an experiment to discover an 
EDM does not depend upon the error in the computed 
enhancement factor being small. Estimates of the size of 
limit that may be set by a francium fountain experiment 
operated at existing accelerator facilities may be found 
in Appendix A of Ref. [27]. 

A key part of any EDM measurement is controlling any 
systematic error that, like the effect of the EDM, reverses 
sign when the electric field is reversed. Such effects arise 
naturally in EDM experiments because in the rest frame 
of the atom, due to the atom's nonzero velocity in the 
applied electric field, the atom sees (SI units) a magnetic 
field 

B mot = \vxE, (1.5) 

which changes sign when E changes sign. The interaction 
of the atom's magnetic moment with that motional mag- 
netic field (Zceman effect) then also changes sign when 
the direction of E is reversed. The rotation of the atom 
by this motional magnetic field, necessarily perpendicular 
to the electric field, cannot by itself generate a rotation 
about the electric-field axis. A priori, one would thus 



TABLE III: Calculations of the enhancement factor R for the 
ground state of heavy alkali atoms. 



Alkali 


R 


Year 


Reference 




24 


1966 


[29] 


Rb 


24.6 


1985 


[30] 




25.68 


1994 


[31] 




25.74 


2008 


[32] 




119 


1966 


[29] 




114.9 


1985 


[30, 33] 


Cs 


114 


1990 


[34, 35] 




120.54 


2008 


[32] 




124 


2009 


[36] 




1150." 


1966 


[29] 


Fr 


910 


1999 


[37] 




894.93 


2009 


[38] 



"This early value did not include a correction for the shielding 
factor of the atomic core. The addition of a shielding correction 
lowers the enhancement factor of all other alkali atoms. 



assume that the motional magnetic field docs not mimic 
an EDM. However, in any practical apparatus, there will 
inevitably be present trace static magnetic fields which 
the atom will explore as it moves; the combinations of 
rotations about these fields and about the motional field 
are of concern. 

In this paper, we present a complete categorization 
of all systematic errors that result from an atom in 
an atomic fountain being subjected to a constant elec- 
tric field, the motional magnetic field, and trace mag- 
netic fields which are static but may vary arbitrarily in 
all three spatial directions. The unitary operator that 
gives the time evolution of any superposition of mag- 
netic sublevels | F, A/) is a complex matrix of dimen- 
sion (2F + 1) x (2F +1), where F is the total spin of 
the hyperfine level. In an atomic fountain, the maxi- 
mum atomic velocity ~ 4 m/s is two orders of magnitude 
smaller than in previous atomic beam experiments [15], 
so the motional magnetic field also is reduced by two or- 
ders of magnitude. In an atomic fountain, atoms pass 
through each point in space twice, once rising and once 
falling; the resulting atom- by-atom reversal of the atomic 
velocity cancels some motional-field dependent systemat- 
ics without a need for a second, antipropagating atomic 
beam. Let us briefly comment on the effect of a possible 
slight nonuniformity of the electric field, in which case the 
atomic trajectory may become more complicated because 
atoms in the level Si/2,F = I + 1/2 are attracted into 
regions of strong electric field by the Stark shift. How- 
ever, in our differential setup, we plan to measure the 
signal twice, with the sign of the electric field inverted, 
and in this case, even a more complicated trajectory is 
conserved and independent of the state of the atom, in 
the approximation that otp 3> in Eq. (1.9). 

Moreover, under the conditions of an atomic fountain, 
the shift of the magnetic sublevel | F, M) due to the ten- 
sor Stark effect is now much greater than its shift due to 
the Zeeman effect. In a constant electric field, the unitary 
operator that gives the time evolution of any state can be 
computed as a scries in reciprocal powers of a dimension- 
less parameter n ~ 100, which represents (roughly) the 
ratio of the shifts. The (l///)-series therefore is rapidly 
convergent. We present this series as the sum of analytic 
integrals over the magnetic fields, times matrices of con- 
stants, for any of the total spins F = 3, 4, and 5, and so 
for the spins of all the alkali atoms of experimental inter- 
est. We place special emphasis on the case F = 4. A key 
result is that a properly constructed observable sensitive 
to an EDM will contain errors that are only even powers 
of l//i. Because contributions of order 1/fi 4 prove negli- 
gible in practice, any experiment only has to control the 
few terms of order 1 //1 2 . 

The coordinate system used to describe the atomic 
fountain is shown in Fig. 1. The atom remains on a 
single vertical axis, the y axis. The direction of the elec- 
tric field is horizontal and defines the z axis. Therefore, 
the motional magnetic field is parallel to the x axis. The 
effective Hamiltonian within a manifold of states of total 
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spin F is the sum of several pieces. The contribution of 
the electron EDM is 



H 



EDM = R ( — d e ) F ■ E , 



(1.6) 



where F is the total angular momentum of the atom, g? c 
is the electron EDM, and R is the enhancement factor. 
The contribution of the Zccman effect is 



Hz = Ms 9f F ■ B , 



(1.7) 



where (Ib = h\e\/ (2m e c) is the Bohr magneton and g F 
is the Lande g factor for the manifold. The Stark effect 
contributes to each level | F, M) the energy shift 



-(a + a FM )E 2 



(1.8) 



where ao is the scalar polarizability (which is indepen- 
dent of F and M) and a F M is the tensor polarizability. 
The tensor polarizability splits [39, 40] into the sum 



O-fm = a F 



3 M 2 — F(F - 
F(2F - 1] 



1) rp 



(1.9) 



where all dependence on the magnetic quantum num- 
ber M is explicit. The parts of ao and of olfm that 
are independent of F z contribute to a global shift of the 
whole hypcrfine manifold and thus introduce no change 
in the atomic state other than a global phase. We may 
therefore drop them in solving for the time evolution of 
the atomic state. In doing so, we are well aware of the 
fact that, while the shifts do not affect the state, they still 
have a large effect on the motion of the atom. In general, 
the global shift of the hypcrfine manifold implies that 
an atom accelerates as it enters an electric field. Fur- 
thermore, a parallel beam of atoms defocuses as it enters 
the electric field, because atoms are pulled towards the 
high-field region at the edges of the plates. Such defocus- 
ing can be compensated by a suitable set of electrostatic 
lenses [41-43]. 

Isolating the terms that are relevant for the quantum 
dynamics (mixing within the F -manifold), this gives an 
effective Hamiltonian 



Hs — A s E F 2 



where As is the constant 



As = - 



3a£ 



2F(2F- 1) 



(1.10) 



(1.11) 



The total effective Hamiltonian H t then is the sum of the 
tensor Stark term, the EDM term, and the Zccman term, 



H t = A s F 2 E 2 +[i B g F F-B-Rd c F z E z 



(1.12) 



Contributions to the effective Hamiltonian from the mix- 
ing of different hypcrfine levels or from terms in the Stark 
effect of order E z (the hypcrpolarizability [44]) are neg- 
ligible. A remark about the units in the Hamiltonian is 



necessary here. We have deliberately scaled the quanti- 
ties in the Hamiltonians so that the total angular mo- 
mentum F as well as the projection quantum number F z 
of a quantum state are dimcnsionlcss. To this end, we 
have absorbed % as the unit of the angular momentum 
into the respective moments and into the constant As, 
using Eq. (1.3) for the EDM and using the usual defini- 
tion of [IB- 

The atoms enter the electric field and spend a time T 
within it. In this paper, we do not consider effects due 
to a continuous rise of the electric field from zero, but 
model the rise as a step-function. In principle, the fi- 
nite size of the transition region over which the electric 
field rises from zero to its full value will contribute to the 
time evolution operator in a rapidly convergent series in 
powers of the small time atoms spend in the transition 
region. The resulting systematics can be controlled by 
reducing static magnetic fields only in the transition re- 
gion, whose vertical length will only be a few centimeters, 
without having to reduce the static magnetic fields every- 
where in the whole fountain, whose vertical length will 
be about a meter. In order to simplify the calculations, 
we also introduce a dimensionless variable t that runs 
from — t to t, with r = 1/2, while the usual time runs 
from zero to T. We then define a standard electric field 
strength Eg by 



TA S 



(1.13) 



in order to render the Schrddingcr equation dimension- 
less. (Note: to use this equation the units of As have 
to be J/(V/m) 2 .) The atom now evolves according to 
H(t)ip = idip/dt, where the Hamiltonian H(t) is ob- 
tained from H t given in Eq. (1.12) by an appropriate 
scaling, 



H(t) = el (t) F 2 + {3(t)-F + <j F e z (t) F z 



(1.14) 



Now H(t) is dimensionless, and the (dimensionless) co- 
efficients are 



E z (t) 
E s 



g F [Ib B(t) 
A S E 2 



a F 



d c R 



FA S E S 
(1.15) 

The time-dependent Hamiltonian (1.14) is the basis of 
the entire derivation that follows. When the electric field 
is a constant in time, we define a time-independent pa- 
rameter 



(1.16) 



which for a typical experimental setup has values in the 
range fi ~ 100. 
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II. TIME-ORDERED PERTURBATION 
THEORY 

A. Hamiltonian and time evolution 

Specializing the above general statements to an atom 
with a defined quantum number F, we now study the 
effective Hamiltonian in Eq. (1.14) for cesium ( 133 Cs), 
whose upper hyperfine level has F ~ 4. As a simple 
model, consider the fountain of Fig. 1 with the atoms 
confined to the y axis, with an electric field which might 
vary in magnitude but is always parallel to the z axis, 
and a magnetic field which in the atom's rest frame varies 
only slowly in time. We define t as a dimensionless time 
variable. The ordinary (dimensional) time i, as measured 
by a clock, is related to t as follows, 

.-(i-i). 

The turning point of the atoms is at t = | T. The atoms 
enter the fountain at t = 0, which corresponds to the 
scaled time variable t = — \, and they leave the fountain 
at t = +5. In the following, we employ the more general 
notation 

-T<t<r, T =\- ( 2 - 2 ) 

This use of the symbol r allows us to recognize time- 
symmetries in the calculation. For example, the motional 
magnetic field in the rest frame of an atom is odd around 
t = because the velocity changes sign there, while lab- 
oratory static magnetic fields are even in t. 

The Hamiltonian (1.14) gives rise to a time evolution 
operator. Complete knowledge of this operator would al- 
low us to propagate an initial state through the fountain 
and obtain the final state, which determines the observ- 
ables. All effects which are either larger or about the size 
of the target sensitivity of 2 • 10~ 50 C m have to be cal- 
culated if they could possibly mimic an electron EDM. 
At this stage, it is useful to recall that the traditional 
unit for the electron EDM is e cm; the conversion to SI 
units is lCm = 1.602 ■ 10 _21 ecm. At the target sensi- 
tivity, the corresponding dimensionless parameter op in 



our Hamiltonian (1.14) has the value up = 4 x 10 9 for 
cesium. 

In the model calculation we consequently need to find 
all effects that would lead to an EDM-like signal of a 
magnitude greater than that of up. Without the EDM 
term, the Hamiltonian reads as 

H (t)=el(t)F*+P(t)-F. (2.3) 

The time-evolution governed by this Hamiltonian is char- 
acterized by the equation 

H (t)U (t)=i^U (t), (2.4) 

with the initial condition U(—t) = 1. 

We first split this time evolution operator U(t) into 
a diagonal part V(t) and a remainder term W{t). 
For W(t), we define an equation of motion and 
write it as a time-ordered exponential W(t) = 
Texp(-i J l T h(t')dt'), where h(t) is defined in Eq. (2.11) 
below, and T denotes the time ordering (recall the sym- 
bol T denotes instead the total time the atom spends 
in the fountain). Because W(t) is expressed in terms 
of highly oscillatory integrals, for a constant magnetic 
field it is amenable to an expansion in the parameter 1 j \x 
where /j, is defined in Eq. (1.16), with coefficients that 
depend on time-integrals over various components of the 
magnetic fields. 

In the absence of magnetic fields in the x and y direc- 
tion, the Hamiltonian Ho(t) becomes 

H^t)=e 2 z (t)F^+p z (t)F z . (2.5) 

This Hamiltonian is diagonal in the hyperfine manifold, 
and its time-ordered exponential 

V{t) = Tcxp (-\ j H^dA (2.6) 

is a diagonal (2F+1) x (2F+1) matrix acting on a 2F+1 
dimensional subspace of M levels inside the manifold of 
a given total spin F , which can be written as 



V(t) = diag 



/ 

exp 
cxp 



-iJ(F(5 z (t') + F\l(t'))dt' 

— T 

-i}{{F-l)P z {t') + {F-lYel(t>))dt> 



exp 



-if {{-F)p z (t') + {-Ffel(t')) dt> 



(2.7) 



/ 



where diag denotes the diagonal matrix generated by the given matrix elements. It is now our aim to solve the 
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full problem by writing the full time-evolution operator Uo(t) as the product of the time-evolution operator for the 
diagonal part V(t) and the remainder term W(t), i.e. 

U (t) = V(t)W(t). (2.8) 

The required equation of motion for W(t) can be found by using this product for U(t) in Eq. (2.4), which gives 



H Q (t) V{t)W{t) = i 1 (V(t)W(t)) , H (t)V(t)W(t) - i^W) = i V(t)^^- 



Multiplication with V 1 (t) from the left leads to the equation 
. dW(t) 



at 



= h(t)W(i) , h{t) = V- l (t)H (t)V(t) - iV-\t) 



dV(t) 

at 



(2.9) 



(2.10) 



The matrix representation of this new Hamiltonian h{t) for general F as a (2F + 1) x (2F + 1) matrix can be written 
in terms of the following, somewhat self-explanatory notation, where the 2F + 1 elements on the diagonal are written 
in the middle, the 2F entries of the sub-diagonal on the left, and the 2F entries of the super-diagonal on the right, 

/ d F f(t)exp\-i(2F-l) f_ T el(t')dt'] d F f* (t) exp \+i(2F - 1) f_ T 4(f) 



h(t) = diag 



d*._i/(t) exp \-i(2F - 3) f_ T e 2 z (t') dt'] d F ^f*(t) exp f+i(2F - 3) f_ T e 2 z (t') dt' 



d F f(t) exp [+i(2F - 1) /* r e 2 z (t>) dt'\ d F f*(t) exp \-i(2F - 1) f_ r e 2 z (t') dt' 





(2.11) 



Here, f(t) is a single complex function that contains all 
information about the magnetic fields in the atom's rest 
frame. 



P z {t')dt') . (2.12) 



f(t) = (p x (t)+ip y (t)) exp 



The coefficients d F need to be explained. They are de- 
fined as the corresponding entries (rational numbers and 
square roots of rational numbers) of the matrix represen- 
tation of the x component of the total angular momen- 
tum F, which is F x = (F + + F_)/2. These entries are 
symmetric with respect to the transformation M — > —M, 
which is why we have been able to reduce the notation 
to the quantities dy[ with M = 1, . . . , F on the sub- and 
super-diagonals. 

A central point of our proposed scheme is to require 
that the integral over the square of the cumulative scaled 
electric field strength seen by the atom along its path 
should be equal to an integer multiple k e of n, 



l{t) dt = k ( _ 



(2.13) 



The adjustment to an integer value of k e has been used in 
a prototype experiment [45, 46]. According to Eq. (2.7), 
if the evolution of the hyperfine sublevels were exclusively 
given by the electric term (no motional or stray magnetic 
fields), then for fc e even, the atomic wave function would 
return to its initial quantum state if the quantization 
condition (2.13) is fulfilled. 



By measuring the temporally constant stray magnetic 
field seen by the atom, and by purposely applying a small 
additional magnetic field over the electric field plates, it 
is possible to also enforce a condition for the integral of 
the z component of the magnetic field seen by the atom 
in the experiment, 



P z {t) dt = kf}n. 



(2.14) 



Here, it is the assumed that hp is tuned to be integer, and 
we note that hp = is preferred; this value corresponds 
to small overall magnetic fields. Assuming that these 
adjustments can be achieved perfectly, the diagonal time- 
evolution operator V(t) at the time t = t when the atom 
leaves the fountain takes the simple form 



V(t) = diag 



(_l)(^-i) 2 fc e 



(_l)i^ ^ 



V(-i) ( - 



FYk e 



(-1) 



-Fk R 



(2.15) 



/ 



which for even k e + kp reduces to the unit matrix. 

We express the remainder W(t) of the time evolu- 
tion operator as a time-ordered exponential in terms of 
a Dyson series. As we are only interested in its effect at 
the time t = t, when the atoms leave the fountain, we 
write 



W(t) = 1 + 



OO 

E 



(2.16) 
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The terms have the form 



r ti 



tn-1 



w n = (-i)" / J ■■ j h{tl) h{t 2 ). . .h{t n ) dt n ... dt 2 dh . 

— T —T —T 

(2.17) 

For the finite time interval t G [— r, r] considered in this 
work, the series for W(t) is in fact convergent. 

The basic idea of our theoretical analysis is as follows. 
The terms w n describe the deviation of the quantum dy- 
namics from the idealized result (2.15) due to the mag- 
netic field (2.12). Analytic estimates of the terms w n 
will allow us to gauge the magnitude and significance 
of individual EDM-mimicking signals with respect to our 
target sensitivity. An asymptotic expansion of the terms, 
based upon a separation of integrands into fast oscillating 
electric-field terms and slowly oscillating magnetic-field 
terms, leads to those estimates. 

The asymptotic expansion is best illustrated by con- 
sidering the first order term 



w\ = (-i) 



ft(ii)d*i 



(2.18) 



The nonzero entries of the (2F + 1) x (2F + 1) matrix w\ 
are of the form 



J = 



/(i) exp in / e 2 z {t')dt')dt 



(2.19) 



where we use n to represent the integers in the exponents 
in Eq. (2.11), and f(t) is from Eq. (2.12). We envisage 
a situation where |e z (i)| 2 is on the order of a hundred. 
The phase in the above integral is therefore highly oscil- 
latory and the integral can be expressed in terms of an 
asymptotic expansion, as it is for example explained in 
Ref. [47]. This expansion is obtained by repeated inte- 
gration by parts of Eq. (2.19) as 



J 



>/-^(*')dt' 



-fit) (2.20) 



3 in/* r eJ(t')dt' 



/"(*) " 2/'(t) 



This equation expresses the time integral over an arbi- 
trary integrand function f(t) as a function of surface 
terms which need to be evaluated at the upper and lower 
limits of the time interval within the fountain. In our 
case, f(t) is given by Eq. (2.12) and describes the de- 
pendence on the magnetic field seen by the atom as a 
function of the scaled time t. The expansion allows us 
to integrate out the fast oscillations due to the applied 
strong electric field and to concentrate, in the formula- 
tion of the quantum dynamics, on the dependence due to 
the residual magnetic fields (static and motional). 

We can further simplify the expression for the expan- 
sion by noting that the exponential 



exp 



\{t')dt' 



(2.21) 



evaluated at t = — r, assumes the value of unity. In all 
calculations where we employ this expansion, the integer 
n in the exponentials in Eq. (2.11), which has one of the 
values 2F - 1, 2F - 3, . . ., -(2F - 1), is odd. With the 
adjustment of the integral in Eq. (2.13), the exponential 
can be evaluated at t = r to be 



exp 



\{t')dt' 



(-!)» 



(2.22) 



If in addition we assume the electric field to be constant, 
the formula for the asymptotic expansion of the integral J 
then is written as 



J = 



fit) exp 



— T 

i 



in / e 2 z (t')dt' 



dt 



„((-U*-/(r)-/(-r)) 
+ ^ ? ((-l)' Se /V)-/ / (-r)) 



+ 



n 3 [i 3 



((-l) fe =/"(r)-r(-T))+0( M - 4 ) 



(2.23) 



Such an expansion provides us with a tool to analyze 
the relative magnitude of the relevant physical effects. 
In our units used in this work the parameter /i for ce- 
sium is about fi rj 120, and for the envisaged sensitiv- 
ity of an electron EDM the dimensionless EDM strength 
in Eq. (1.14) is a F = 4 x 10~ 9 . Therefore, all EDM- 
mimicking effects of order 1/u 4 4 • 10 -9 have to be 
considered and eliminated in order to reach the target 
accuracy. In practice the actual coefficients of terms of 
order 1//Z 4 are small enough that their contribution is 
negligible, and for a suitable observable all effects of or- 
der l//i" for n any odd integer can be proved to can- 
cel, leaving only two terms of order to compute 
and to control. Using asymptotic expansions like that of 
Eq. (2.23) for the integrals appearing in W(t), we are 
able to express the time-evolution operator 

U q (t) - V(t) W(t) = Texp f -i / H Q (t')dA (2.24) 



when the electric field is constant in terms of an asymp- 
totic series in inverse powers of fi, which has the structure 



Mr) 



1 



(2.25) 




f(n) {t) rW(r) 



dt 



where the % are complex (2F + 1) x (2F + 1) dimen- 
sional coefficient functions that can be expressed in terms 
of integrals of the nth derivatives of the stray magnetic 
field function f(t) or its complex conjugate f*(t), and 
in terms of powers thereof. This functional dependence 
is schematically indicated in Eq. (2.25) using the curly 
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brackets and the multi-index {n,m, k,£}. We have car- 
ried out the calculation to fourth order in l//z, for F = 3, 
F = 4, and F = 5. The individual expressions are too 
long to include in this paper, but the calculation is en- 
tirely based on the asymptotic expansion technique de- 
scribed herein. 

The expansion of Eq. (2.25) has applicability beyond 
the case when the electric field is constant. Provided 
e z (t) is approximately constant on [— r, r] there exists 
a transformation of the time variable that produces an 
equivalent problem where the electric field is exactly con- 
stant and where the unitary transformation is exactly of 
the form of Eq. (2.25), with a perturbed function /. As 
far as the cancellation of systematic effects in an electron 
EDM experiment, nothing has changed; only the numer- 
ical values of surviving systematics are perturbed. 

Before we use this expansion to evaluate the specific 
terms of the time-evolution operator, we want to take 
the time and discuss the observable to be used in the 
experiment [45, 46]. The choice of the observable also 
identifies the EDM-mimicking effects. 



B. Defining the observable 

Without magnetic fields in the x and y direction, the 
only effect of the movement through the fountain for the 
atoms would be a rotation of the quantum states by a 
complex phase [see Eqs. (2.7) and (2.15)]. This follows 
because the phase adjustment of the magnetic field inte- 
gral as well as the integral of the square of the electric 
field relevant to the Stark shift, sets the complex phase 
equal to a multiple of it. For k t + hp even, the complex 
phase is a multiple of 27r and so the atoms are precisely 
rotated back into the initial state when they leave the 
fountain. The difference between the Hamiltonian with 
the EDM term given in Eq. (1.14) and the Hamiltonian 
without the EDM term given in Eq. (2.3) is that the 
presence of an EDM leads to a small additional rotation 
around the z axis, and so the complex phase would no 
longer be zero (or equal to an integer multiple of 2tt). 
Therefore we need an observable that is sensitive to ro- 
tations about the z axis. 

In the case of 133 Cs, the valence electron is in the 6S' 1 / 2 
state, so the atom can be in either of the two hypcrfinc 
levels F = 3 or F = 4. The atoms are prepared in a 
state [4, M). After the atoms leave the fountain, the 
number of remaining atoms in the F = 4 hypcrfine level 
is measured and compared the total number of atoms. 
Some atoms can also transition into the F = 3 hypcrfinc 
level of the 651/2 state. 

Let us denote the state in the F = 1+1/2 hyper- 
fine level in which the atom is prepared before it enters 
the electric field as |$o)- After the atom exits the elec- 
tric field, it is analyzed by optical pumping with a laser 
propagating in the z direction and tuned to a P state 
with spin F' = F, where the laser is linearly polar- 
ized. We are interrogating quantum transitions with re- 



spect to quantum states whose quantization x and y axes 
are tilted by a rotation an angle 8 with respect to the 
z axis. The states on which we are projecting thus are 
the states R z (8) | F, M) , where R z (8) is a suitable rota- 
tion operator. In this paper, we use the notation R u (a) 
to indicate an operator that rotates a state on which it 
acts (active representation) about the axis denoted by 
u and by an angle a that is positive if the rotation is 
in a positive sense about the axis u as determined by 
the usual right-hand rule. Under these assumptions, we 
have R z (0)\F,M) = exp(-iM0) \F,M). The probabil- 
ity that the atom will be found in the dark state of the 
upper hyperfine level is our basic signal and is given by 

F 

S{0)= Yl Pm\(FM\R z (6)U(t)\$ )\ 2 . (2.26) 

M=-F 

where the constants pu can be chosen to be any of the 
sets in Tables I or Table II. Here, U(t) is the time- 
evolution operator of the full Hamiltonian (1.14), 

U(t) =Texp f-ij H(t)dtj . (2.27) 

A small rotation of the atomic states by an EDM can 
then be detected by taking two measurements at ±8 and 
forming the difference. Thus, the observable P{8) is given 
by the equation 

P{6) = S{6) - S{-8) 

F 

= Y,Pm(\(FM\R z (+8)U(t)\$ )\ 2 (2 28) 

M=-F 

-\(FM\R z (-9)U(t)\$ )\ 2 ) , 

where because we are using a linearly polarized laser the 
probabilities pm have the symmetry pm = P-m- This 
observable is measured again with the direction of the 
electric field reversed. The resulting difference is only 
sensitive to effects that, like that of an EDM, are odd 
under the reversal of the electric field direction. 

The currently most promising alkali atoms for an ac- 
tual experiment are cesium and francium. For cesium, 
the natural isotope 133 Cs is suited best. It has nuclear 
spin I = 7/2, and so the ground state 66*1/2 has the 
two hyperfine levels F = 3 and F = 4. The proposed de- 
tection scheme requires to use the energetically higher 
lying state F = 4. For francium the enhancement fac- 
tor R is about nine times larger than in cesium, mak- 
ing this atom particularly attractive. The francium iso- 
tope 221 Fr, which has I = 5/2, can be obtained from 
actinium sources and the dynamics of the F = 3 hyper- 
fine level would be investigated there. Due to a higher 
nuclear spin / = 9/2, leading to a higher total angular 
momentum F = 5 to be used for the dynamics, the fran- 
cium isotope 211 Fr, is currently the most promising atom 
to study. It can be obtained from accelerator sources, 
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such as ISOLDE and TRIUMF, and its half-life of 3 min- 
utes is sufficiently long for a practical experiment. Yields 
of francium at existing and planned facilities (CERN 
ISOLDE and TRIUMF) and a discussion of the transfer 
mechanism from the beam to the magneto-optical trap 
can be found in Refs. [48-53]. 

In the calculation of P(9), it is possible to perform a 
basis transformation which leads to a number of simpli- 
fications. We define as our new basis states the following 
superpositions of the sublevels ±M, which are obtained 
from the original basis states \FM) by a rotation about 
an angle tt /2 around the y axis, 



1 

71 
i 

71 



R y {n/2) (\FM) + \F — M)) , (2.29) 
R y (v/2)(\FM)-\F -M» . (2.30) 



Here R y (Tr/2) is the corresponding rotation operator, and 
we index the new states by their dependence on M and 
suppress the subscript F. These new states may also 
be used for the projections carried out in the calculation 
of P{6). 

The new states have useful symmetries under two 
transformations. Let us define the matrix B, which has 
the entries 



Ba = (-1) 



i+l 



(2.31) 



where the indices run from 1 to 2F + 1. This matrix is 
trivially generalized to F = 3, 4, and 5: the B matrix 
is a diagonal (2F + 1) x (2F + 1) matrix where the first 
matrix elements on the diagonal run +1, — 1, +1, , . . . The 
states \sm) are eigenstates of B with eigenvalue +1, while 
the states |om) are eigenstates of B with eigenvalue —1. 

The new states also have a symmetry under exchange 
of magnetic quantum numbers from M to —M. This 
exchange can be described by transformation by a second 
(2F + 1) x (2F + 1) matrix S that has entries 



Sij — 5 it (2F+l)-j 



(2.32) 



Both \sm) and \qm) are eigenstates of S with eigenval- 
ues (-1) F+A/ . 

In terms of these new states, and in the case where 
Pm = P-Mi we can rewrite P(9) as 



P(6) = J2p M (\(s M \R z (6)U(T)\$ )f 

\( Sm \rI(o)u(t)\^ )\ 

J2PM(\(aM\R z (d)U(T)\<S> )t 
\(a M \Rl(0)U(T)\$ o )\ 



M=0 



(2.33) 



M=l 



where we have used the fact that R z (—6) = R\(8). 

The great utility of the inversion symmetry becomes 
clear in the next step, as it allows us to express the dif- 
ference of the measurements at ±9 instead in a simpler, 
product form [see Eq. (2.35) below]. We define the "even" 
part of a matrix A as 

0+{A) = ±(A + SAS), 

which is even under conjugation with S, 

SO+{A)S = ±S(A + SAS)S = 0+{A) , 

because SS = 1. The "odd" part 

0_{A) = ±(A-SAS) 

is odd under conjugation with S, 

SO-(A)S= \S{A- SAS)S = -0-(A). 

Due to the identity [we use R z = R z (0)] 



|<V> 2 |^A|Vi)| 2 - \{ih \RtA\ ^)| 2 = -4Re(Vi \0_{A^)RM{^ \R\0 + {A)\ V>i> , 



(2.34) 



we have the following expression for P(6), 



P(0) 



4 Re 



4 Re 



Pm <$o |0-(t^(r))| Rt(0)s M ) (R z (0)sm \0+{U(r))\ <K ) 

Vf=0 
F 

E Pm <*o \0-(UHr))\ R\{6)a M ) {R z {9)a M \0 + (U(r))\ $ ] 



M=l 



(2.35) 



This is an improvement on Eq. (2.33) in that we have expressed the difference between the measurements at ±9 as the 
real part of a product: of the part of the time-evolution operator evaluated at t = r that is even under conjugation 
with S\ and a part which is odd under conjugation with S. In addition, we have separated the projection onto 
intermediate states that are even and odd under conjugation with B, namely |sjw) and Iom), respectively. 
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For k e + kp even, we can use W(r) instead of U(t) because V(r) is the identity matrix. For k e + kp odd, the 
effect of V(t) is a minus sign for each of the matrix elements containing a^j for F = 4, which as we will see later do 
not contribute to an EDM or EDM-mimicking signal. We are therefore safe to just use W(t) instead U(t). Using 
Eq. (2.16), P(0) can now be rewritten as 



P(6) 



= Re 



M=0 
F 



Pm($0 O-l 



00 \ t 



Rt{9)s M )(R z (6)s M 



M=l 



Ri(6)a M )(Rz(0)aM 



0+(l + ^w„ 



(2.36) 



Here, use is made of the fact that 0_(1) = 0. The observable now being known, we can employ the asymptotic 
expansion to determine the matrices 0±(w n ) which describe the time-evolution of the atoms through the fountain to 
identify EDM-mimicking effects. We begin by determining the terms in the first order of the Dyson series. 



C. First term of the expansion of W 



We continue to use the Hamiltonian h(t), which does not contain a possible EDM, in order to isolate all effects 
that could mimic the presence of an EDM. We obtain the unitary transformation accurate to order l//i 2 , for which 
as we shall see it is necessary to consider some terms at fourth order in time-ordered perturbation theory. Knowing 
the terms in the unitary transformation to order l//x proves to be sufficient to describe all EDM- mimicking effects 
in a suitable observable with error of order l//z 4 . The calculational scheme for our observable P(9) implies that we 
have to determine the 5-even part of the matrix w\, denoted as 0+{w\), and the 5*-odd part, denoted as 0-{w\), 
separately. We present these characteristic calculations in some detail and start by giving the matrix for Oj r {w\), 
which for general F has the form 



0+(t«i) = diag 



id F jRc[f(t)]cxp -i(2F-l)Je 2 z (t')dt 



— T 
I 



id F _! / Rc[/(t)] cxp -i(2F-3) / e 2 z (t')dt' 



-id F /Re[/(i)]exp +i(2F-l)fe 2 z (t')dt' 



dt 
dt 

dt 




-ieZ F fRe[/(i)] exp +i(2F-l)je 2 z (t')dt' 

— r — r 

T . t . 

-id F -i/Re[/(f)] exp +i(2F-3) / e 2 z (t')dt' 



idf /Re[/(t)] exp -i(2F-l)Je 2 z {t')dt' 



dt 
dt 

dt 



J 



(2.37) 

where the (If are the corresponding entries in the matrix F x . A simple calculation reveals that Sw\S can be obtained 
out of W\ by exchanging f(t) and f*(t). The entries of 0+(wi) are thus f(t) + f*(t), which is twice the real part 
of /(*)■ 

Similarly, in 0_, the difference of the function f(t) from Eq. (2.12) and its complex conjugate gives the imaginary 
part of f(t). There is an additional minus sign for the super diagonal because of the different order of / and /*. As 
a result of the calculation, 0-(wi) is found to be 



( d F ] Im[/(t)] cxp \-i{2F-l)Jtl{t')dt' 



0-(wi) = diag 



dt 



d F J Im[/(t)]exp +i(2F-l)/ e 2 z (t')dt' 



d F _i/Im [/(<)] exp -i(2F-3)fe 2 z (t')dt' dt -d F -if Im[/(t)] exp +i(2F-3)/ e 2 z {t')dt' 



d F J Im[/(t)]exp +i(2F-l)/ e 2 z (t')dt' 



dt 


I 



dp J Im[/(t)] exp -i(2F- 1)/ e 2 z (t')dt' 



dt 

/ 

(2.38) 



If n = 



— (2F — 1), . . . , 2F — 1 denotes the prefactor in the exponentials of the matrix elements of 0+(uii), then 
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when t z is constant wc find by asymptotic expansion with 
error 0(/i~ 3 ) we have 

[O+Mn = ~ [Mf(-r)} Rc[/(r)](-l) fe '] 

+i§^[Re[/'(-r)]-Re[/'(r)](-l)^] . 

(2.39) 

The function / is given by 

f(t) = [p x (t)+ip y (t)]e- iB V, (2.40) 

with 



B(t) = / 0,(0 dt'. 



(2.41) 



We can simplify results by using the time-symmetry of 
fields in an atomic fountain. Because the atom is assumed 
to rise and fall along the y axis, it passes through the 
same static magnetic fields when both rising and falling. 
In the rest frame of the atom the magnetic fields applied 
this way are time-even, i.e., even around the time t = 0. 
Similarly the electric field applied in atom's rest frame, 
due to the atoms motion through a static electric field 
that we have assumed is parallel to the z axis, is time- 
even. The motional magnetic field applied in the atom's 
rest frame, which arises due to the Lorentz transforma- 
tion of the electric field and in our geometry is always par- 
allel to the x axis, is time-odd around t — because the 
atomic velocity changes sign when the atom falls while 
the direction of the electric field does not. 

We can split the magnetic field along the x axis, (3 x (t), 
into its static and time-even part x e (t) and its motional 
magnetic and time-odd part x (t). Thus we have the 
following facts: 

• x Q (t) = fi x (motion; t) is time-odd, 

• x e (t) = /3 X (static; t) is time-even, 



y e (t) = P y (static; t) is time-even 
e z {t) is time-even around t = 0, 



• z e (t) = P z (static; t) is time-even, 
and in addition B(t) = (3 z (t') dt' = kpir. 

Again specializing to the case of a constant electric 
field e z , to order l//z we find 

Re[/(-r)] - Re[/(r)](-l) fe « = Re [&(-t) + i^(-r)] 

(2.42) 



Separating the magnetic field in the x direction into 
its static and motional part and employing the time- 



symmetries in the ideal fountain, we obtain the result 

&(-t) - (-l) fc * +fc ^x(r) 

= x e (-r) + x (-t) - (-l) k <+ k e [ Xe ( T ) + Xo (t)] 



2x (~t) for k e + k/3 even , 
2x e (— t) for k e + kp odd . 



(2.43) 



To get the contribution of order 1 //i 2 , we need the deriva- 
tive 



= ^ [&(*)+ i&(t)] e- iB ^ 



(2.44) 



= [(3' x (t)+W]e- iB W 

which allows us to obtain 

Re[/'(-r)] - Rc[/'(r)](-l) fe = = y e (-r)z e (-r) 

-(-l) k ^y e (-r)z e (-r) 

+ <(-t) + x' a (-r) (-l) k ^x' e (r)x' (r) . 

(2.45) 

Because the time derivative an odd function of time is 
even, and the derivative of an even function is odd, we 
get 

y e (-r)z e (-r) (-l) k ^y e {-r)z e (-T) 

+ x' e (-r) + x' (-r) (-l) k ^x' e (r)x' (r) (2.46) 

{2x' e (— r) for fc e + kp even , 

2y e (-r)z e {-T) + 2x' (-t) for k e + k p odd . 

We also define a standard form for the expansion of the 
matrices 0+(w n ), which takes the form 

k=k ^ j 

where the \ are matrices of dimension (2F + 1) x 
(2F + 1) whose entries are real, field-independent con- 
stants that depend on F. The expressions are also 
real and can have a rather complicated dependence on 
the magnetic field /3(t), as well as integrals, derivatives 
and of powers thereof, but are independent of F. The 
sum over j is introduced to effect a natural separation 
of the terms of a given order in n and k according to 

specific symmetry properties of the field-dependent func- 

(i) 

tions as explained below. 

In first order in n, where the subscript n denotes the 
expansion order for the time-ordered perturbation the- 
ory, only a single term with j = 1 is required, and we 
can suppress the superscript j in our notation; the sum 
over j becomes necessary only for the second and higher 
orders (n — 2, 3, . . . ). Low-order terms in n will nonethe- 
less contribute to high orders in l//i. Corresponding to 
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each n is however a lowest-order nonvanishing contribu- 
tion in l//i. While one's initial guess would be that the 
sum over the inverse powers of 1/fj, starts at k = n, in 
fact this sum starts at the value ko defined by 



ko = k {n) 



n/2 for n even, 
(n+ l)/2 for n odd, 



(2.48) 



or in an alternative notation as the ceiling of n/2 (i.e., the 
smallest integer larger than or equal to n/2). This more 
complicated rule arises from the cancellations among the 
various exponential factors in various integrands, as ex- 
plained in more detail below. 

For 0-{w n ), we define an analogous series by 



o_K)= J2 



k—kn 



E 





(2.49) 



where we now use M. to denote the respective matrices 
of real field-independent constants; B for the respective 
magnetic-field dependent, real expressions; and where ko 
is as defined in Eq. (2.48). The matrices M. and M as well 
as the expressions Q and B can be shown to be real to all 
orders. Because the only imaginary factors in Eqs. (2.47) 
and (2.49) are the explicit powers of i, it is straightfor- 
ward later to take the real part of the observable P{9). 
Suppressing for simplicity of notation the magnetic 

(i) 

field dependence of G^ k , the leading contributions 
to 0+{w\) can be written 



o+( Wl ) = --jv 1a (f) g 1A 



(2.50) 



with the functions 



2x (— t) for k e + kp even, 
2x e (— t) for k t + kp odd, 



(2.51) 



and 



n,2 - 



2x' e (—T) for k e +kp even, 

2y e (-T)z e (-T)+2x' (-T) for fc.+ fc^odd. 

(2.52) 
4 have 



The coefficient matrices J\f for cesium with F 
the form 



M,i(4) 





/ -a/2/7 


V2/7 \ 




-Vu/w o 


Vu/10 




-a/2/2 


V2/2 




-VE o 


V5 


diag 


a/5 


-V5 




V2/2 


-V2/2 




Vli/io 


-y/U/10 




V2/7 


-V2/7 




V 


J 



(2.53) 



and 



M,2(4) = -diag 



/\/2/49 \/2/49\ 

/14/50 VT4/50 

a/2/6 \/2/6 

vt \/S 

a/5 \/5 

a/2~/6 v^/6 

/14/50 \Ai/50 

a/2/49 a/2/49 




(2.54) 



In the notation introduced in Eq. (2.49), we find for 
0-(w) 



0-(w 1 ) = —M 1 , l (F)B 1A 

+ \m 1 , 2 (F) B 1 . 2 + 0( f i- 3 ), 
The magnetic field dependence is given by 



(2.55) 



Si, 







for fc e + &8 even, 



2y e (—r) for fc c + fc^ odd. 



(2.56) 



and 



_ |2j/^(-r)-2a;o(— r)^e(-r) for fc e +fc,3 even, 
1 — 2x e (— r)z e (— r) for k e +k/3 odd. 

(2.57) 

For cesium with F = 4, the coefficient matrices are given 



A^i,i( 4 ) = - dia 8' 



-a/2/7 





-a/2/7 \ 


-vTi/10 





-Vu/10 


-\/2/2 





-a/2/2 


-V5 





—a/5 


a/5 





a/5 


\/2/2 





a/2/2 


\/l4/10 





a/14/10 


a/2/7 





a/2/7 










(2.58) 



and 



-Mi, 2 (4) = diag 



/ -a/2/49 





a/2/49 \ 


-a/14/50 





a/14/50 


-a/2/6 





a/2/6 


-V5 





a/5 


— V^5 





a/5 


-a/2/6 





a/2/6 


-a/14/50 





a/14/50 


-a/2/49 





a/2/49 


V 





/ 



(2.59) 



For atoms with F = 3 and .F = 5, the coefficient matrices 
are given in Appendices E and F of Ref. [27] , respectively. 
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III. HIGHER-ORDER CALCULATION 



For the diagonal part of order 1/fJ,, we have 



A. Second term in the expansion of W 

In the second order of time-dependent perturbation 
theory, we need to investigate the matrix W2 , which takes 
the more complicated form 



Wo 



h(ti)h{t 2 ) dt 2 dii 



(3-1) 



— T J —T 



and where there are more nonzero entries in the resulting 
matrix. Due to the presence of two factors in the inte- 
grand, each with a phase factor, one's initial guess would 
be that the resulting matrices are at least of order 1 //i 2 . 
However, terms of lower order in l//i appear because of 
cancellations in the exponentials. For example, consider 
the terms on the diagonal of 0+(w2)- The structure of 
these integrals is 



J 



jn[E{ti)-E(t 2 )] 



x [/(i 2 )/*(ti) + /(ti)/* (t 2 )] 



(3.2) 



dt 2 d<i , 



where n = 2F-1, 2F-3, . . . , -(2F- 1) is an odd integer 
and we have used 



E(t) 



4(f) dt' 



(3.3) 



to shorten the notation. When the electric field is con- 
stant, we have E(t) = [i 2 (t + r); and carrying out the 
asymptotic expansion for the first integral leads to 

^=S/j /(ti)|ad * i+ ^£^ l/(ti)|adti 

+ -^( 2 |/(-r)| 2 - /( T )/*(- T )(-l)*« 

(3.4) 

whose lowest contribution is seen to be of 0(1/ ^l), 
not 0(l//i 2 ). In general the contribution of lowest or- 
der in l//i in the expansion of w n is not given by the 
order in perturbation theory but rather only by the ceil- 
ing of n/2; hence the introduction of kg in Eq. (2.48) 
and in the asymptotic series expansions of Eqs. (2.47) 
and (2.49). 

Bringing the expressions for 0+(w 2 ) into the standard 
form of Eq. (2.47) yields 



+ (™ 2 ) = -i^4?(F)Sg 



+ ^E^(f)^] + o(m- 3 ). 



(3.5) 



< 1 i(4) = -diag 



/ -2/7 \ 
-29/70 
-4/5 
-7/2 

10 
-7/2 
-4/5 
-29/70 
V -2/7 / 



and for the off-diagonal part we have 

F(F + 1) 



(3-6) 



M%(F) 





F(F + 1) Q 



V 



(3-7) 



7 



where we use the dots to indicate that all other entries of 
the (2F+1) x (2F + 1) matrix N 2 [2 }(F) arc zero. For the 
diagonal part, the magnetic field dependent coefficients 
are 



»2.1 



(3-8) 



{xl(t) + x 2 (t) + y 2 e (t)} dt, 
and for the off-diagonal part 

gj$ = 2(-l) k > fj[x 2 e (t) + * 2 o(t) - yl(t)]™sg(t) 

+ 2x (t)y e (t)sing(t)\dt, (3.9) 



where we have defined the new function 

ri 



9(t) 



P z (t')dt'. 



(3.10) 



o 



which is odd in time. The diagonal term of order l//i 2 is 



/ -2/49 \ 
-443/2450 
-16/25 
-11/2 
-10 



A/- 2 ( ^(4)=diag 



with 



»2,2 



-11/2 

-16/25 
-443/2450 
-2/49 / 



for k e + fcg even, 



(3.11) 



2x 2 (-r) 

2x\(-t) + 2j/ 2 (-t) for k t + k p odd. 



(3.12) 



16 



There is a second term, with superscript j 
coefficient matrix 

/ 



< 2 2 } (4) 



2\/7/35 
V7/5 
2\/7/35 VlO 

V7/5 -10 

vTo vTo 

-10 y/7/5 

VTo ooo 2V7/35 

v/7/5 
\ 2\/7/35 / 

(3.13) 

Here and in the following we refrain from using the "diag" 
notation for matrices whose structure is beyond tridiag- 
onal. Blank entries are understood to denote zeros. The 
field-dependent terms contain the magnetic fields squared 
and so are even when the electric field is reversed, 



2, with a has the form 



\ 



- 

F(F + 1) Q 



F(F+1) 



\ 



■J 



(3.18) 

where we use the notation established in Eq. (3.7). The 
associated magnetic field dependence is electric-field even 
and given by 



£fi = 2(-l)W x e {t)y e (t) cos g(t)dt. (3.19) 



fl(2) _ 

y^2 — 



for k e + hp even, 



x\{-t) - vI(—t) for k e + hp odd. 



(3.14) 



A third term with j 
given by 



3 is odd in the electric field and 



<2 ) (4)= 





n/7/105 

^720 

\/l0/2 



-\/7/105 
-V7/20 







-VT0/2 




T0/2 
-V7/20 

-\/7/105 



/TO/ 2 
V7/20 
%/7/105 
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with 



(3) 
»2,2 



a; e (-T)x (-T) . 



/ 

(3.15) 
(3.16) 



For 0_(w 2 ) we find, written in the standard form (2.49), 



0.(w 2 ) = ±Y i M$(F)B$ 



(3.17) 



Interestingly, 0_(w2) has no term of order l//x on the 
diagonal. In the terms on the diagonal, the exponential 
factor is the same as the exponential factor of Eq. (3.2), 
but with the plus sign changed into a minus sign; when 
the asymptotic expansion is carried out, the contribu- 
tions of order l/fi now cancel instead of add. Off the 
diagonal there are still terms of order 1/fJ,. These con- 
nect only states with \M\ = 1 and the coefficient matrix 



There also is an electric-field odd contribution with the 
same coefficient matrix, i.e. 



m?Uf) = m£Uf), 



(3.20) 



(2) 

and with B\ [ given by 



82} = -2(-l)^y x (t)x e {t)smg(t)dt. (3.21) 



There are four matrices that occur at the order l//i 2 . We 
have 



/ -2/49 \ 
-243/2450 
9/25 
9/2 

, (3.22) 



M 2 %A) = diag 



9/2 
9/25 
243/2450 
2/49 / 



where the magnetic field dependence is a mix of terms 
that arc electric-field even and clcctric-ficld odd, 



°2,2 



z e (t) [x 2 e (t)+x 2 (t)+y 2 e (t)] dt 

-2 f_ T x' Q {t)y e (t) dt for k t + kp even, 
+2 f_ r x (t)y' e (t) dt for k e + kp odd, 

(3.23) 
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and 



B. Calculation of the observable 



/ v^/105 

V7/20 
s/7/Wh vTO/2 

V^/20 

y/W/2 



-vTO/2 



-VW/2 
-y/7/20 
-vT/105 
-\/7/20 
— n/7/105 o o 

(3.24) 

where the magnetic-field dependence is electric-field odd, 
only: 



V 



/ 



B 



(2) 
2.2 



x (-r)y e (-r) . 



(3.25) 



The next matrix is 



/ -4-/F/35 

-2-v/7/5 



M5(4)= 



4^/35 
2\/7/5 
2v / 10 



-2VT0 


ooo -2V10 

-2V7/5 
2VlO -4VT/35 
2V7/5 
4V7/35 







with the electric-field even dependence 



(3) 
2,2 



for fc e + fc^ even, 

x e (—T)y e (—T) for fc e + k/3 odd; 



(3.26) 



(3.27) 



and the last matrix, which only acts on states with \M\ 
1, is given by 



(■ 



Mi%F) 



\ 



F(F + 1) 

-F(F+1) 



where the dependence is electric-field even, 



a 



J for fc e + A:^ even , 

1 x e {-r)y e (-T) for fc £ + fc^ odd . 



(3.28) 



(3.29) 



This concludes the contribution from w 2 . The terms of 
order l//x 2 that derive from third and fourth order in 
perturbation theory are given in Appendices C, D, G 
and H of Ref. [27], respectively. 



At this stage, we can keep all terms of order 1 //x 2 in 
our observable P{6) by truncating Eq. (2.36) to 



-P(6)/4 



Re 



^ Pm(®0 \0^(w 1 +W 2 + W 3 + W 4 )' I \ Rz 



sm) 



M=0 



x (R z s M \0 + (l + w 1 + w 2 + w 3 +w 4 )\ $ ) 



^ Pm($o \0-(wi+ w 2 +w 3 + w 4 )^\ R z 



M=l 



x (R z a M \O + {l+w 1 +w 2 +w 3 + w 4 )\<S> ) 

(3.30) 

Simplifications of this expression are possible using the 
two symmetry operators S and B. Through these sym- 
metries, the states can be classified into orthogonal sub- 
spaces. Matrix elements between states from different 
subspaccs then vanish due to the orthogonality. This, 
however, requires us to specify the initial state and the 
total angular momentum. Otherwise, the behavior with 
respect to the S and B symmetries cannot be determined. 

Consequently, we now specialize to F = 4 and to one 
of the initial states 

|$ ) = ^(|4,4) + |4 J -4)), (3.31) 

|$„> = ;^(|4,2> + |4,-2». (3.32) 

Both of these have the same behavior with respect to the 
two symmetries B and S, i.e. 



B|$o> = (+l)|*o> 
5|$o> = (+l)|*o> 



(3.33) 
(3.34) 



Let us examine the matrix element in the first line of 
Eq. (3.30) 

(^$0, \0-(lVi + ^2^ + ^3^ + w 4,)^\RzSm) ■ 
b=+l b=-l b=+l 6=-l b=+l b= + l 

and examine how the pieces transform under B. We see 
that the state 0-(wiY\R z sm) belongs to the eigenspace 
of B with eigenvalue b = —1, while |$o) has eigen- 
value +1, and so the overlap of these states vanishes. 
Similarly the overlap involving 0-{w 3 )^ vanishes, and so 
the matrix element as a whole reduces to 

-> (<S> \O-(w 2 +w^\R z s M ) ■ (3.35) 
An analogous analysis for the other matrix elements al- 
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lows to reduce P(6) to 



-P(0)/4 = Re 



Y1p m \0-(w 2 +w 4 ) i \R l 



M=0 



x (i? z s M IO+(l+w;2+W4)|4>o> 



+ »( $ o (O-^i+wa) 1 ! R z a M ) 



M=\ 



x ( J R 2 a M |0 + (wi+W3)|*o) 



(3.36) 



We now consider the symmetry of the various contribu- 
tions under S. The real and the imaginary part of the ro- 
tation matrix R z transform differently, i.e. SRc(R z )S = 
+Re(R z ) and SIm(R z )S = -Im(R z ). Recalling that 
each of the initial states |$o) under consideration have 
eigenvalue +1 under S, we can eliminate further terms 
from P{ff) and find 



P(9)/4 



Re 



-i^ p M (^o\0-{w 2 +w 4 ) ii \lm(R z )sM) 

M=0,2A 

x (Rc(R z )s M \0+(l +w 2 + io 4 )| $ ) 
+i Pm ( $ o \0- (w 2 + w 4 ) f I Re(R z )s M ) 

M=l,3 

x (lm(R z )s M |0+(1 + u; 2 + w 4 )| $0) 
-i Pm{$0 \0-(wi +w 3 y\lm(R z )a M ) 



M=2,4 



x (Re(i? z )a A/ |0 + (wi + tu 3 )| *o) (3-37) 
+i pm($o \0-(wi +w 3 )i\Rc(R z )a M ) 



M=l,3 



x (Im(i? z )a M |0+(wi + 103) I $0) 



The last step in the simplification is to take the real part. 
Because the total expression carries a pre-factor of i, only 
those products of matrix elements contribute that con- 
tain an additional imaginary unit. From Eq. (2.47), wc 
recall that the term in order \j\i m of 0+(w n ) is imagi- 
nary as \ n+m is imaginary, and from Eq. (2.49) that the 
term in order l//j, m of 0-(w n ) is imaginary as i™+ m+1 
is imaginary. So, in a product of these two matrices the 
total phase is \ n ++ n -+ m ++ m -+ 1 _ Here the subscripts + 
and — refer to a contribution from O+ and 0_, respec- 
tively, and n refers to the order of a term in perturbation 
theory and m refers to the power in at which the 
contribution contributes. In each of the products of ma- 
trix elements in Eq. (3.37), the phase i"++ n - is always 
real because n+ and n_ are either both odd or both even. 
The remaining phase [ m ++ m -+ 1 [ s imaginary only when 
the order in l//i, equal to m++m_, is even. Thus, in the 
final result for P{9) only even powers of 1/fx survive. In 
fact this result is perfectly general: the same arguments 



applied to the full expression of Eq. (2.36) show that the 
asymptotic expansion for P{6) contains only even powers 
of 1/V, for all powers however high. 

Therefore though we have kept only terms of order 
1/fj, and I//1 2 in the unitary transformation in writing 
Eq. (3.30), and therefore while one might expect the er- 
ror in the resulting expression for P(9) to be of order 
l//i 3 , the error is in fact of order l//x 4 . We confirm 
this result by comparing the result for the observable in 
Eq. (3.43) below with a numerical solution for the prob- 
lem, as shown in Fig. 3, where the difference between 
the calculation here and a numerical solution multiplied 
by /j, 4 is plotted. The difference converges to a constant 
for fi ^S> 1. Similar tests have been performed to ensure 
that the expansion of the time-evolution operator U(t) 
includes all effects up to including order I//1 2 . 

//A [fj 4 10" 43 C m] 

A 



4.75 



4.55 



20 



40 



60 



71 



FIG. 3: Plot as a function of /j, of the scaled difference 
A = fJ, 4 (pA ~ P°Anum); here T° A is defined by Eq. (3.47), and 
the difference is taken between the result from the asymptotic 
expansion based on Eq. (3.43) and the result from a numer- 
ical calculation PXnum- The scaled difference approaches a 
constant as fi — > 00 and thus confirms the error estimate. 

At this point we can go back to our full Hamiltonian 
including an EDM, which as we recall from Eq. (1.14) 
has the form 

H{t) = e 2 z (t) F 2 Z + /3(i) ■ F + ap e z (t) F z . (3.38) 

The electron EDM d e enters the dimensionless EDM cou- 
pling of = —d c R/(F As Eg); it is assumed to be very 
small and we only need to include this effect in first or- 
der. We use a time-evolution operator for the EDM part 
of the Hamiltonian in Eq. (1.14), 

W = l-iJ a F e z (t)F z dt = l + iDF Zl (3.39) 

where we have defined 

D = -a F J e z (t)dt, (3.40) 

where o\f has been defined in Eq. (1.15). The EDM term 
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adds the D term to our observable. 



different initial states |$o) evaluates to 



P(0)/4 
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In this expression, we now take a closer look at the terms 

(2) 

and realize that the important difference between B 2 2 
and i?2 3 2 i s tnat ^22 reverses sign when the direction of 

the externally applied electric field is reversed while B 2 \ 
stays the same. The signal in which we are interested 
is D, which also changes sign when the electric field is 
reversed. So, if we measure P(9) twice, once with the 
electric field along the z direction, and again with the 
electric field in the — z direction, and subtract the results, 
then all electric field even terms are canceled in P(9). For 
the electric field odd part of P(6), we find 



V°(9) = -{P(9) +E - P(9)_ E ) 

= 4^p M (lm{R z )a M \^o) {®o\DF s 
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Here, the electric-field odd part B 2 2 
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The set of probabilities pm with M = 0, 1, 2, 3, and 4 

can be chosen from among the sets in Table I by choosing 

the frequency of the analysis laser. The angle 9 can be 

varied continuously by changing the inclination of the 

axis of linear polarization of the analysis laser. 

In the expressions above the angle 9 is in the range 
(2) 

[0,7r/2]. The term in B 2 2 is even in 9 about the mid- 
point 9 = 7r/4, and the term in B 2 X 2 ^ and in the electron 
EDM contribution D are both odd. The odd linear com- 
bination for data for different angles 6, 
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(3.47) 



both cancels any contribution from B 2 2 and maximizes 
the sensitivity of what remains to D. The even linear 



combination isolates B^\. Once isolated, B 2 L 2 ' can be 
tuned to be of small magnitude because it depends on 
the vertical component of the static magnetic field at the 
entrance to the electric field plates, which can be tuned 
by varying a current in a nearby coil. 

(2) 

Once B 2 2 has been tuned to be of small magnitude, 
For F = 4, the observable P(0)[$o] corresponding to two canceled, or both, the remaining systematic B^ 1 ^ can 



(1,0) 



M=0,2,4 

^[M^B^f+M^B^l] \lm(R z )a M ) 

+ o( P - 4 ). 
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be canceled by taking the right linear combination of 
data for different initial states, as shown by Eqs. (3.43) 
and (3.44). This systematic can also be isolated, and 
then tuned to be of small magnitude by varying the ver- 
tical component of the static magnetic field, in this case 
at a location away from the entrance to the electric field 
plates. 

A big advantage of the proposed observable is that 
there are no EDM mimicking effects of order 1 //i 3 , which 
had it been present would have been relevant on the pro- 
posed level of sensitivity. It is not yet necessary to cancel 
or control EDM-mimicking effects of order 1/^t 4 , partic- 
ularly since in practice, our experimental accuracy will 
be limited by other systematic effects such as magnetic 
Johnson noise, as for example discussed in Ref. [54]. 



IV. EXPERIMENTAL REALIZATION 

So far, we have considered a rather idealized model 
of an atomic fountain. In order to make sure that an 
actual experiment is not hindered by systematic effects 
arising from more realistic conditions, we now consider 
systematic effects outside the model. 

The effect on V ^ of a linear vertical gradient in the 
component of the magnetic field in the ^-direction is ex- 
plored numerically in Fig. 4. Because of the parabolic 
rise and fall of the atoms while within the electric field, 
such a gradient results in the atomic rest frame as 
a y component of the magnetic field that is of the form 
j3 y (static) (a + bt 2 ) for constants a and b. In Fig. 4, we use 
atomic data for the atom 133 Cs and find that depends 
linearly on both a and b. The numerical calculation which 
produces Fig. 4 has also been verified against the ana- 
lytic result found in Eq. (3.45) for the EDM-mimicking 

effect B^f. 



-4 




-1.0 



FIG. 4: Plot of a numerical result for the signal function given 
in Eq. (3.47), denoted as for 133 Cs, depending on a and 
b. The a and b parameters describe the magnetic field in the 
y direction as f3 y (t) = f3 y (static) (a + bt 2 ). 



A similar numerical analysis is shown in Fig. 5. of the 
dependence of on a deviation of Sk e and Skp from 
their respective phase adjustment to integer multiples fc e 
and kp of 7r. The results suggest that the deviations 
should not exceed ~ 10~ 3 for the error in the signal to 
remain within the limits set forth by our target accuracy 
of - 10~ 9 . 




FIG. 5: Plot of a numerical result for the signal of Eq. (3.47), 
for 133 Cs depending on the deviations 5k t of k e and Skp 
of kp from the phase adjustment in Eqs. (2.13) and (2.14), 
respectively. 

In a real as opposed to an ideal fountain, the common 
point of state preparation and analysis does not occur at 
the edge of a step-function rise of the electric field. In 
practice both will occur at a point where the electric field 
is essentially zero, well below the entrance to the electric 
field plates where over a transition region whose verti- 
cal extent is the order of the plate spacing the electric 
field ramps from zero to a constant value. Between the 
common point and the beginning of the transition region, 
the atoms will be in the presence of a stray magnetic field 
that will effect a small rotation of the atomic state. In 
the transition region the combination of a significant mo- 
tional magnetic field and a stray magnetic field, acting 
while the tensor Stark splittings are still small, generates 
among other effects an additional rotation about the z 
axis that can mimic the rotation effected by an electron 
EDM. All such effects influence the atom both on its way 
into the full electric field and on its way out. 

The additional effects can be studied by defining time- 
evolution operators for each of these processes. For the 
transition region, we may expand the time-evolution op- 
erator in terms of the short (dimcnsionless) transition 
time H, which is small against the overall time T the 
atom spends in the fountain. The rotation that occurs 
(about some arbitrary axis) before the atom enters the 
transition region may be expanded in the small rotation 
angle S. We end up with three parameters l//i, K, and S, 
which in the natural units of this work are all dimension- 
less, small, and in a practical apparatus roughly of the 
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same order of magnitude, ~ 1/fj,. We can thus system- 
atically model EDM-mimicking effects up to a desired 
order in the small parameters. Terms of third order are 
relevant on the proposed level of sensitivity. A detailed 
discussion of this expansion, which follows a formalism 
akin to the calculation of the time-evolution operators 
considered here, remains beyond the scope of the current 
paper. Technical details, which will be given elsewhere, 
reveal that the dominant terms in H and in <5 do not 
mimic an EDM and can thus be discarded. The remain- 
ing effects can be satisfactorily described on the level of 
our proposed target accuracy. 

One might also speculate about the effect of a depar- 
ture of the atomic motion within the electric field from a 
straight line, as well as the effects of misalignments of the 
axes of propagation or of polarization of the laser used 
for state preparation and analysis. The first of these ef- 
fects can be included within our approach by adding an 
electric-field odd but time-even motional magnetic field 
in the y direction. We have carried out a calculation 
of this effect whose detailed discussion again is beyond 
the scope of the current paper. We find that only the 
field-dependent functions change but not the coefficient 
matrices. Our proposed cancellation mechanism is based 
solely on the coefficient matrices and so not affected by 
a change in the values of the field-dependent functions. 
Misalignments of the direction of propagation of a laser 
or of its axis of polarization can be modeled by repre- 
senting the actual laser as being an ideal laser rotated by 
a small angle; this small angle can be treated as a small 
parameter that affects the time evolution of the atoms 
just as do l//x, H, and 5. We have not found any effects 



that could possibly represent insurmountable difficulties 
at the level of our proposed target accuracy; the design 
of the experimental apparatus is in progress. 



V. FRANCIUM AND CESIUM SYSTEMS 

The availability of accelerator sources has made possi- 
ble the use francium (half-life 3 minutes) instead of ce- 
sium in the atomic fountain. The main advantage of 
francium is that the enhancement factor R is about 9 
times larger in francium than in cesium (see Table III). 
Furthermore, certain francium isotopes, especially 211 Fr, 
have a higher total angular momentum, namely F = 5, 
which reduces the relative size of the EDM-mimicking 
effects. Moreover, the higher tensor polarizability also 
leads to \x being larger and thereby speeds the conver- 
gence of the formalism presented here. 

Here, we present the results of our time-ordered cal- 
culation for F — 3, which covers the case of 221 Fr, and 
for F = 5, which covers 211 Fr. Only the coefficient ma- 
trices jV and M. differ from those for F = 4; the corre- 
sponding field-dependent functions B and Q are identical. 
In the asymptotic expansion of the unitary transforma- 
tion U to order 1//J 2 , those matrices not given for gen- 
eral F in Sec. II are gathered in the Appendices E, F, 
G, and H of Ref. [27]. Here, we only present the re- 
sults for our proposed observable T' (9)[^o\ for four ini- 
tial states |<I>o) for the two isotopes; results for some 
other initial states are given in Appendix I of Ref. [27]. 
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The variable D has been defined in Eq. (3.40) and incor- 



porates the effect of the EDM in view of its dependence 
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on ap defined in Eq. (1.15). Equations (5.1) and (5.2) 
show that when F is odd, initial states consisting of su- 
perpositions of states with equal \M\ do not exhibit the 
symmetry in 9 around tt/2 that when F is even could 

(2) 

be used to cancel the systematic effect B\^- Measure- 
ments on the states \F, 0) are shown by equations (5.3) 

and (5.4) to be sensitive only to only to B\ which can 
therefore be measured and tuned to be small by changing 
the vertical component of the static magnetic field at the 
entrance to the electric field plates. Any residual effect 
can be separated from the values of D and of the other 
systematic jE?2 2 either by combining data for different 
angles 9, or by altering the laser frequency and combining 
data for different sets of values pm- The sets accessible 
for F = 3 and for F = 5 may be found in Appendix B of 
Ref. [27]. 

VI. CONCLUSIONS 

We set up the formalism for the theoretical descrip- 
tion of the quantum dynamics of an alkali atom within 
an atomic fountain designed for an EDM experiment. 
The low velocity of the atoms inside a fountain reduces 
the motional magnetic field, which arises as the Lorentz 
transformation of the applied electric field, by a fac- 
tor of 100 compared to experiments on thermal atomic 
beams. In an atomic fountain the quantization axis is 
defined by the direction of the externally applied strong 
electric field, and is not defined by a magnetic field; we 
are therefore free to greatly reduce all magnetic fields 
and their attendant systematic errors by use of extensive 
magnetic shielding and nulling coils. Compared to many 
previous experiments these two features suppress effects 
that mimic the presence of an EDM, because such ef- 
fects scale linearly with both the motional magnetic field 
and any stray field. We describe a theoretical calculation 
that identifies the remaining EDM-mimicking effects and 
devise schemes to eliminate them. 

A crucial part of our formalism is writing the time evo- 
lution operator from time-ordered perturbation theory in 
terms of an analytic expansion in the inverse number of 
electric-field induced Rabi oscillations within the hyper- 
fine manifold. When the magnetic fields in the x and y 
directions are zero, the time-evolution operator for a sin- 
gle hyperfine manifold is a diagonal matrix of phases V, 
which reduces to a simple rotation of the atomic sys- 
tem provided the cumulative electric and magnetic field 
phases of Eq. (2.13) and Eq. (2.14) are tuned to integer 
multiples of n. 

The perturbing effects of stray magnetic fields in the x 
and y directions as well as the motional magnetic field 
in the x direction, still without an EDM term, can be 
treated when the electric field is constant by expanding 
the correction term W(t) in the formula for the time- 
evolution operator U(t) = V(t) W(t) in an asymptotic 



expansion in inverse powers of \i = e 2 , where /z/(2-7r) de- 
scribes the number of Rabi oscillations in the hyperfine 
manifold induced in the absence of perturbing fields. Un- 
der typical experimental conditions a value of /x on the 
order of one hundred can be obtained and so expansion 
is therefore rapidly converging. The expansion of W to 
order 1/fi? requires consideration of terms of fourth order 
of perturbation theory; we present the expansion in terms 
of constant matrices multiplied by analytic functions of 
the stray and motional magnetic fields. 

We then define an observable P(9), which is a func- 
tion of the angle of inclination of the linear polarization 
of the laser used to analyze the final state of the atom. 
This observable is sensitive to a rotation of the initial 
state about the electric field axis, as would be induced 
by the presence of an electron EDM. A transformation 
of the final state into pieces symmetric and antisymmet- 
ric with respect to two transformations B and S defined 
in Eqs. (2.31) and (2.32) proves that many contributions 
cancel. Critically it is proven that in the asymptotic ex- 
pansion of P only even powers of l//x appear, so knowl- 
edge of the time development operator out to order l//i 2 
suffices to compute P with error l//i 4 , not l//^ 3 . 

After taking the difference for opposite signs of the 
electric field, besides the effect of an electron EDM only 
two systematic errors survive, which can both be isolated 
and canceled by combining data for polarization angles, 
laser frequencies, and initial states; moreover both sys- 
tematic errors once measured can be tuned to be of small 
magnitude by imposing additional small magnetic fields. 
The systematics start intrinsically smaller than in previ- 
ous experiments because of the smaller velocity of atoms 
in an atomic fountain, which reduces the magnitude of 
the motional magnetic field, and because in an atomic 
fountain with a large electric field a magnetic field is not 
required to define a quantization axis, so stray compo- 
nents of the magnetic field can be suppressed to the limit 
provided by the surrounding magnetic shielding. 

The proposed apparatus is expected to detect an EDM 
on a level of 2 x 10 -50 C m, or better. The limit is due to 
higher-order effects and due magnetic Johnson noise [54] 
from the materials used in the apparatus. 
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